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ABSTRACT 

The key properties of the wave propagation theory in the magntosphere of radio pulsars 
based on the Kravtsov-Orlov equation are presented. It is shown that for radio pulsars 

with known circular polarization and the swing of the linear polarization position 
angle one can determine which mode, ordinary or extraordinary one, forms mainly 
the mean profile of the radio emission. The comparison of the observational data with 
the theory predictions demonstrates their good agreement. 
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1 INTRODUCTION 

There are three main propagation effects in the magneto- 
sphere of radio pulsars. They arc refraction, cyclotron ab- 
sorption, and the limiting polarization effects. The limiting 
polarization effect is related to the escape of radio emission 
from the region of dense plasma, where the propagation is 
well described within the geometrical optics approximation 
(in this case, the polarization ellipse is defined by the orien- 
tation of the external magnetic field in the picture plane), 
into the region of rarefied plasma, where the polarization of 
the wave remains almost constant along the ray. This pro- 
cess has been well-studied (Zhelcznyakov, f 996; Kravtsov & 
Orlov, 1980) and applied successfully, e.g., to the problem 
of solar radio emission (Zheleznyakov, 1970). 

However, in the theory of the pulsar radio emission such 
a problem has not been finally solved. Above the papers, 
where the position of the transition level rcsc between the 
domains of the geometrical optics and vacuum approxima- 
tions was estimated (see, e.g., Cheng & Ruderman 1979; 
Barnard, 1986), one can actually note only the papers by 
Pctrova & Lyubarskii (2000a) (whore the problem in the in- 
finite magnetic field was considered), and by Wang, Lai & 
Han (2010) (see Petrova, 2003, 2006 as well). But in all these 
papers the equations describing the evolution of the electric 
field of the waves were analysed. This approach cannot make 
the direct predictions concerning the polarization properties 
of the outgoing radiation. 

We use another approach describing the propagation 
of electromagnetic waves in weakly inhomogeneous media, 
i.e., the method of the Kravtsov-Orlov equation that is well 
known in plasma physics and crystal-optics. It gives us the 
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opportunity to write directly the equations on observable 
quantities, i.e., the position angle and the degree of circular 
polarization (CP). Also, the cyclotron absorbtion that nat- 
uraly occurs in finite magnetic field, and the possible linear 
transformation of waves can be easily included into consid- 
eration. 

The preliminary results based on the simple model of 
the magnetic field structure and energy distribution of parti- 
cles flowing in the magnetosphere were already published by 
Andrianov & Bcskin (2010). It was shown that for radio pul- 
sars with known circular polarization and the swing of the 
linear polarization position angle one can determine which 
mode, ordinary or extraordinary one, forms mainly the mean 
profile of the radio emission. Later, the arbitrary non-dipole 
magnetic field configuration, the drift motion of plasma par- 
ticles, and their realistic energy distribution function were 
taken into account as well. The detailed quantitative analy- 
sis of these effects will be presented in our separated paper. 
The goal of this Letter is in qualitative comparison of the 
main predictions of the theory with observational data. In 
our opinion, they are in a very good agreement. 



2 THEORETICAL PREDICTIONS 

2.1 On the number of outgoing waves 

Starting from the pioneer paper by Barnard and Arons 
(1986), people discussed three waves propagating outward 
in the ijulsar magnetosphere. As shown in Fig. 1, for small 
enough angles 6 between the wave vector k and the exter- 
nal magnetic field B two of them, rii and n2, arc transverse 
waves, and the third mode ns corresponds to plasma wave. 
The point is that in the most of papers (see e.g. Melrose 
& Gedalin, 1999) the waves properties were considered in 
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Figure 1. Dependence of the refractive indexes n on angle 9 for 
Ap 2> 1- Tlie lower brancli corresponds to the O-mode 



the comoving reference frame in which the plasma waves 
propagating outward and backward are identical. But in the 
laboratory reference frame (in which the plasma moves with 
the velocity v m c) the latter wave is to propagate outward 
as well. Thus, in reality we have four waves propagating 
outward. 

Moreover, as was demonstrated by Beskin, Gurevich & 
Istomin (1993), it is the fourth wave that is to be consid- 
ered as the O-mode in the pulsar magnetosphere. Indeed, 
as shown in Fig. 1, for dense enough plasma in the radio 
generation domain for which Ap ^ 1 where 
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< 7 >, 
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it is this wave that propagates as transverse one at large 
angles 9 between k and B, i.e., at large distances from the 
neutron star. Here ujp = (47re^nc/mo)^''^ is the plasma fre- 
quency, ric is the concentration of particles, and < 7 > is the 
mean Lorentz-factor of the outflowing plasma. Two waves, 
for which the refractive index n > 1, cannot escape from the 
magnetosphere as at large distances they propagate along 
the magnetic field lines (and due to Landau damping, see 
Barnard & Arons, 1986). 

In the hydrodynamical limit one can easily obtain the 
dispersion curves shown in Fig. 1 from the well-known dis- 
persion equation for the infinite magnetic field (see, e.g., 
Petrova & Lyubarskii, 2000) 
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For 9 <s:9, and for 9:^ 9,, where 
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there are two transverse and two plasma waves, but for 
Ap ^ 1 the nontrivial transformation from longitudinal to 
transverse wave takes place. It means that for Ap ^ 1 the 
mode n4 can be emitted as a plasma wave, but it will escape 
from the magnetosphere as a transverse one. 

2.2 Kravtsov-Orlov equation 

The Kravtsov-Orlov equation describes the evolution of the 
complex angle O = Bi -|- i02, where Oi is a position angle 
of the polarization ellipse, and Q2 determines the circular 



polarization by the relation V = /tanh202. Here I is the 
intensity of the wave, and V is the Stokes parameter. As 
was shown (Andrianov & Beskin, 2010, Beskin & Philippov, 
2011), the Kravtsov-Orlov equation can be rewritten as 
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where the signs corresponds to the regions before/after the 
cyclotron resonance and 
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We would like to note, that in this equations the circular 
polarization is defined as it is common in radioastronomy 
(positive V corresponds to the co - clockwise rotating electric 
field vector for the observer). Here I is a coordinate along the 
ray propagation, and the angle /3(l) defines the orientation 
of the external magnetic field in the picture plane. In the 
geometric optics region (r < rcsc) the O-mode corresponds 
to polarization Qi ^ /3 + S, and X-mode corresponds to 
Oi x^ + S + t:/2. 

Finally, 6^/^/ are the components of plasma dielectric 
tensor in the reference frame where z-axis directs along the 
wave propagation and the external magnetic field lies in xz- 
plane. In comparison with model considered by Andrianov 
& Beskin (2010), equations Q-Q include into considera- 
tion nonzero components Refe^'y/]. It allows us to take into 
account the electric drift of plasma particles. This effect 
just corresponds to aberration considered by Blaskiewicz, 
Cordes & Wasserman (1991). This effect was also consid- 
ered by Petrova & Lyubarskii (2000a), but for the infinite 
magnetic field only. 

Thus, knowing the magnetic field structure and plasma 
properties of the outgoing plasma (i.e., knowing the dielec- 
tric tensor Si'j') one can determine the observable physical 
quantities, namely, the Stokes parameter V — 7tanh202(oo) 
defining the CP and the position angle p. a. = Oi(oo) char- 
acterizing the orientation of the polarization ellipse. This ap- 
proach is valid in the quasi-isotropic case when there are two 
small parameters, i.e., the general WKB condition 1/kL <C 1 
and the condition An/ni,4 <C 1, where An = ni — 714. 



2.3 Main predictions of the propagation theory 

As the refractive index n4 differs from unity, the appropriate 
ordinary mode is to defied from the magnetic axis until 
9 ^ 9f. As was already mentioned, for the O-mode this effect 
takes place until Ap > 1, i.e., for small enough distances 
from the neutron star r < ta, where 
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Here R, P, and B12 are the netron star radius, rotation pe- 
riod (in s), and magnetic field (in 10^^ G), respectively. Ac- 
cordingly, 7100 = 7/100, vghz is the wave frequency in GHz, 
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Figure 2. Propagation of X- and O-modes in the pulsar mag- 
netosphere. The lower curve on each top panel of each picture 
indicates V, on the bottom panel p. a. curve is shown. 



and A4 — A/10^, where A = ric/ncj is the multiplicity of 
the particle creation near magnetic poles (ncj = flB/2nce 
is the Goldreich- Julian concentration). On the other hand, 
the transverse extraordinary wave with the refractive index 
n = 1 (X-mode) is to propagate freely. As the radius ta is 
much smaller than the escape radius rose (Cheng & Ruder- 
man, 1979, Andrianov & Beskin, 2010) 
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one can consider the effects of refraction and limiting polar- 
ization separately. In particular, this implies that one can 
consider the propagation of waves in the region r ~ rcsc as 
rectilinear. 

As a result, as shown in Fig. 2, if two modes are gener- 
ated at the same distance rgcn < from the neutron star 
where j4p(rgon) > 1, the directivity pattern of the O-mode 
is to be wider than that for the X-mode. Hence, it is logical 
to assosiate the X-mode with the core component of the di- 
rectivity pattern. On the other hand, the O-mode is to form 
the conal component. E.g., assuming the constant plasma 
density across the flow (more realistic regime was consid- 
ered by Petrova & Lyubarskii, 2000), Beskin, Gurevich & 
Istomin (1993) obtained for the radio window width Wr 
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For the O-mode, the theory gives two radiation beam 
widths, depending on the distribution of the energy release 
power in height. 

The main new theoretical result based on the analysis 
of Kravtsov-Orlov equations (|4|-([5]l is that for large enough 
derivative d(/3 -I- 5)/dl, the sign of the Stokes parameter V 
is to be determined not by the poorly defined nondiagonal 
components of the dielectric tensor Sx'y' , when Vo/I = ±l/g 
(Zheleznyakov, 1977) and 



but by the sign of the derivative d(/9 + S)/dl. It takes place 
when one can neglect the first term in the r.h.s. of Eqn. Q 
in comparison with the derivative dOi/di ~ d(/3 + S)/dl. 
E.g., forr<resc {\V\ < 1): 
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I 
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Here x — 0,1/ c is dimensionles distance along the ray. 

Hence, for large enough d{j3 + S)/dx ~ 1 (the total turn 
A(/3 + S) ^ 1 within the light cylinder i?L ~ c/Q), and 
for small angle of propagation S <C 1 through the relativis- 
tic plasma (wy/c ~ 1) the Stokes parameter 1/ (141 is indeed 
much larger than Vq/I — ±l/g resulting from standard eval- 
uation. Moreover, numerical calculations (Beskin & Philip- 
pov, 2011) shows that the sign of the derivative d(/3 + 5)/dl 
is opposite to the sign of the derivative dp.a./d<^. 

As a result, one can formulate the following predictions: 

• For the X-mode (Oi ^ f3 + S + tv/2) the theory pre- 
dicts the SAME signs of the circular polarization V and the 
derivative dp.a./d^. 

• For the O-mode (Oi « /3 -I- 5) the signs of the circular 
polarization V and the derivative dp.a./d(j} are to be OP- 
POSITE. 

• For radio pulsar with the tripple mean pulses we have to 
see the O-mode first, than X- and again the O-mode because, 
as was shown in Fig. 2, the O-mode deviates from magnetic 
axis. 

• In general, the trailing part of the main pulse can be 
absorbed (see Dyks, Wright & Demorest, 2010 as well). 

• The leading component can be absorbed only if the 
polarization formes near the light cylinder (rose ~ ^^l). In 
this case the p.a. is to be approximately constant in given 
mode as the magnetic field here is approx. homogeneous 
(A.Spitkovsky, private communication). 

• Statistically, we expect single profiles for the X-mode 
(if only the X-mode is observed) and the double ones for the 
O-mode. 

In the following section we will try to demonstrate that these 
predictions are indeed in good agreement with observational 
data. 



3 COMPARISON WITH OBSERVATIONS 

3.1 Statistics on single- mode pulsars 

In Table [l] we collected more than 70 pulsars from two 
reviews by Weltevrede & Johnston (2008) and Hankins & 
Rankin (2010) for which both the values of p.a. and V were 
well-determined. As was predicted, most pulsars with double 
(D) mean profile corresponds to the O-mode (the opposite 
signs of V and dp.a./d(f>), and the most pulsars with single 
(S) mean profiles - to the X-mode (the same signs of V and 
dp.a./d(j)). Moreover, statistically the O-mode pulsars have 
wider mean pulses than X-mode, their values being in good 



agreement with theoretical predictions ( 10 1-( 12 \ 
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3.2 Pulsars with almost constant p.a. 

In Table[2]we collected ten pulsars with almost constant p.a. 
from the papers by Johnston et al. (2007) ([1] in the refer- 
ence column), Mitra and Rankin (2010) ([2]), and Johnston 
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Table 1. Statistics of pulsars with known circular polarization V 
and position angle swing. Pulsar period P is in seconds, and the 
window width W50 is in degrees. 



Profile 


Os 


Od 


Xs Xd 


Number 


6 


23 


45 6 




6.8± 3.1 


10.7± 4.5 


6.5± 2.9 5.3± 3.0 



Table 2. Pulsars with almost constant p. a. Here P is in seconds, 
P is in 10"^"", and magnetic field B on the surface of neutron star 
is in 10^2 G 



PSR 


P{s) 


P 


B12 


comment reference 


J0543+2329 


0.25 


15.4 


2.0 


+ [3] 


J0738-4042 


0.37 


1.6 


0.8 


+ [1] 


J0837-4135 


0.75 


3.5 


1.6 


[1] 


J1559-4438 


0.26 


1.0 


0.5 


+ [3] 


J1735-0724 


0.42 


1.2 


0.7 


[1] 


B1907-03 


0.50 


2.2 


1.0 


[2] 


J1915+1009 


0.40 


1.5 


2.5 


[1] 


J1937+2544 


0.20 


0.6 


0.4 


+ [1] 



et al. (2008) ([3]). One can note that their rotational periods 
are smaller than 1 s, i.e., the typical pulsar period. Also four 
of them (indicated by plus in the comment column) show the 
"flatting" of the p.a. swing with the frequency decreasing. 
The interpretation of these effects can be easily given bas- 
ing on the estimate of the escape radius rose Q, where the 
polarization of the radio emission forms (Andrianov & Be- 
skin, 2010). As r^Bc/Rh oc P'^l^v'"^^^ , one can see that for 
small pulsar periods and lower frequencies the outgoing po- 
larization forms closer to the Hght cylinder where magnetic 
field of neutron star, as was already stressed, is almost ho- 
mogeneous. In this case the p.a of the outgoing radiation is 
to be approximately constant within the mean pulse. This 
effect was firstly introduced by Barnard (1986), the main 
distinction of our theory is in self-consistent definition of 
rose basing on the solution of equations (|4|-([5|. 

3.3 Pulsars with tripple mean profiles 

Below we analyse the profiles of several tripple pulsars pre- 
sented in review by Johnston et al. (2007) for which both 
p.a. and V are well-determined at any way at one frequency, 
693, 1374, or 3100 MHz. 

For pulsar PSR J0452— 1759 the triple structure is seen 
at 693 MIfz only. For this reason, it is not surprizing that the 
90°-switch of the p.a. from O-mode to X-mode and return 
is seen on this frequency only. The CP is well-determined in 
the trailing subpulse only, where its sign, as was predicted, 
corresponds to O-mode. 

The tripple structure of the pulsar PSR J0738-4042 is 
seen on the p.a. profile at 691 MHz only. Nevertheless, as the 
p.a. is approximately constant at low frequency (see Table 
2), one can assume that the leading part of the profile is 



absorbed. Then, the leading subpulse is to correspond the X- 
mode, and the trailing one - to the O-mode (the appropriate 
negative CP is seen in this subpulse only). 

In pulsar PSR J1559-4438 the tripple structure O-X-O 
is well-seen at all frequencies 691, 1374, and 3100 MHz, but 
in the mean profile the conal component is detected at the 
frequency 3100 MHz only. The negative CP for core subpule 
corresponding to the X-mode is seen at all frequencies. 

The CP of the tripple pulsar PSR J2048-1616, as was 
predicted, corresponds to O-X-O structure. On the other 
hand, the switch of the p.a. is absent, that can be caused, 
in our opinion, by the averaging technique and should be 
checked in the individual profiles as well (one can see the 
infiuence of the averaging technique on the CP profile in 
Karastergiou et al., 2003). 

Thus, as we see, the properties of the pulsars with triple 
mean profiles are in good agreement with the theoretical 
predictions. The central (core) components of the mean pro- 
files, in general, connect with the X-mode propagating freely, 
while the conal parts correspond to the O-mode deflecting 
via the refraction from the magnetic axis. 

3.4 Pulsars with interpulses 

Here the same analysis for pulsars with interpulses from pa- 
per by Keith et al. (2010) is given. The observations were 
made at the frequency 1.4 GHz. 

The main pulse of pulsar PSR J0627+0706 is to be con- 
nected with the O-mode (it has the opposite signs of p.a. 
and V\ The trailing subpulse of the double profile, as is 
well seen from the p.a. swing (the observable subpulse cor- 
responds to the first half of the S-shape curve) is absorbed. 
In the interpulse the circular polarization is not high enough 
to determine modes. 

The central part of the main pulse of pulsar PSR 
J1549— 4848 is formed by the X-mode (the same signs of 
the p.a. and V). In the leading and the trailing parts of the 
main pulse the CP is too low to say anything about the conal 
component. The interpulse has the triple structiure, the se- 
quence of the modes, as was predicted, being O-X-O. This is 
clear not only from the negative-positive-negative sequence 
of the Stokes parameter V but from the singularity of the 
p.a. at the phase = —85°. 

The main pulse of pulsar PSR J1722— 371 is formed 
mainly by the X-mode. The CP of the interpulse is too low 
to determine the modes. 

The double main pulse of the pulsar PSR J1739-2903 
is to be interpreted as the sequence of the O- and X-modes. 
It is clear both from the negative-to-positive change of the 
Stokes parameter V and from the 90° jump in p.a. The trail- 
ing (O-mode) component is absorbed because, as for PSR 
J0627-I-0706, the leading O-mode subpulse corresponds only 
to the first half of the S-shape curve of the p.a.. The inter- 
pulse is formed by the X-mode. 

In the main pulse of pulsar PSR J1828— 1101 the value 
V is not high enough to determine the mode. The interpulse 
is formed mainly by the O-mode, the trailing subpulse being 
absorbed (here again the p.a. of the visible pulse corresponds 
to the first half of the S-shape curve) . 

Thus, for radio pulsars with interpulse the interpreta- 
tion on the ground of the theory under consideration is rea- 
sonable as well. 
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3.5 Milliseconds pulsars 

Finally, we consider the millisecond pulsars which detailed 
polarization characteristics were presented recently by Yan 
et al., (2011). The observations were made at the frequency 
1369 MHz. 

Pulsar PSR J0437— 4715 has the multiple mean profile. 
Its leading and the trailing parts (i.e., the conal component) 
definitely corresponds to the O-mode. On the other hand, 
both the 90° jump in p.a. and the change of the Stokes 
parameter sign in the core component shows that the central 
subpulse consists not only of the X-mode, but of the O-mode 
as well. 

The leading and the trailing components of the double 
profile of pulsar PSR J1022+1001, in agreement with the 
prediction, are formed by the O-mode. The peak on the p.a. 
curve in the center of the mean profile could be connected 
with the core X-mode, but its intensity is too low to be seen 
on the mean pfofile. 

Pulsar PSR J1045— 4509 has the triple main profile. It 
has different signs of the circular polarization in the core and 
conal components. But the position angle swing is irregular, 
which prevent us to determine the modes. 

In pulsar PSR J1600— 3053 the core component is defi- 
nitely formed by the X-mode (the same signs of the Stokes 
parameter V and the derivative Ap.a./A(j)). Then, the conal 
component is to be connected with the O-mode as its p.a. 
curve locates approx. 90° higher that that for the core com- 
ponent. The circular polarization is too low to confirm this 
point. 

The circular polarization of the central part of the dou- 
ble profile of pulsar PSR J1603— 7202 is different from the 
conal ones. The irregular character of the p.a. swing does 
not allow us to determine the mode. 

Pulsar PSR J1643— 1224 has the single main profile 
which can be easily interpreted as standard O-X-O sequence 
with the absorbed trailing part. Indeed, the polarization of 
the leading part definitely corresponds to the O-mode (the 
different signs of V and dp.a./di/)), while the trailing part is 
to be connected with the X-mode (the same signs of V and 
dp.a./Acj)). Moreover, the derivative dp.a./A(j> is much larger 
in the trailing part. This implies that the visible pulse cor- 
responds the first half of the S-shape curve. 

Pulsar PSR J 1713-1-0747 has multicomponent profile 
containing both modes. It has different signs of the circu- 
lar polarization for different branches of the p.a. curves. 

Pulsar PSR J1732— 5049 has actually zero circular po- 
larization which does not allow us to determine the modes. 
Nevertheless, one can see that the p.a. of the core and conal 
parts belongs to the different modes. 

The profile of pulsar PSR J1909-3744 is quite similar 
to PSR J1643— 1224. Hence, it can be easily interpreted as 
O-X-O sequence with the absorbed trailing part as well. 

Radiation of pulsar PSR J2129— 5721 is to be connected 
with the core conponent formed by the X-mode (the same 
signs of V and Ap.a./A(j>). Two subpulses of the mean profile 
can be easily explained by the central passage through the 
directivety pattern (see Fig. 2). 

In addition, the double profile pulsars PSR J0613— 0200, 
J0711— 6830 have the same signs of the Stokes param- 
eter V in both subpulses. But it is impossibe to con- 
nect them with the O-mode because of the irregular 



swing ofthe p.a.. Finally, pulsars J1024-0719, J1730-2304, 
J1744-1134, J1824-2452, J1857-h0943, 2124-3358, and 
J2145— 0750 have irregular structure or actually zero cir- 
cular polarization which does not allow us to analyse their 
properties. 

Thus, in that cases when the p.a. swing and the Stokes 
parameter V are regular enough, the main properties of the 
mean profiles can be easily interpreted within the theory 
under consideration. 
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